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Abstract 

A recently developed treatment of stochastic processes leads to the construction of a potential 
landscape for the dynamical evolution of complex systems. Since the existence of a potential 
function in generic settings has been frequently questioned in literature, here we study several 
related theoretical issues that lie at core of the construction. We show that the novel treatment, 
via a transformation, is closely related to the symplectic structure that is central in many branches 
of theoretical physics. Using this insight, we demonstrate an invariant under the transformation. 
We further explicitly demonstrate, in one-dimensional case, the contradistinction among the new 
treatment to those of Ito and Stratonovich, as well as others. Our results strongly suggest that 
the method from statistical physics can be useful in studying stochastic, complex systems in general. 
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I. INTRODUCTION 



Stochastic differential equations (SDEs) were first introduced in studying Brownian mo- 
tions originally discovered in biology l| 



2J. Other developments in complex systems, from 



Darwin's evolutionary theory [3(, the dynamics of gene regulatory networks [4|, the cycle 



;uiatory 
igm y, 



kinetics in biochemical networks [5( , to the landscape paradigm |6|] , all have suggested that 
SDE is a useful and effective mathematical description for a wide range of dynamical pro- 
cesses {2, 7|. A distinct difference between the SDEs from biology and those widely studied 
in classical physics, usually called Langevin equations, is the absence of the detailed balance 
condition in the former. Hence the existence of a potential landscape with both local and 
global meanings has been frequently called into doubt [8|. During the study of the robust- 
ness and stability of a genetic switch involved in phage lambda, a bacterium killing virus, 
one of us has discovered a mathematical structure embedded in the SDE. This structure al- 
lows a direct quantification of the global stability of systems with fluctuations. Such unique 
structure leads a transformation of a SDE to a mathematical form familiar to theoretical 
physicists. The transformed SDE has a potential function, a friction matrix, a transverse 
force, and a random noise, together with a relationship between the noise and the friction 
similar to the fluctuation-dissipation relation [9], n]] . These four dynamical components can 
be obtained constructively from a SDE. 

The mathematical structure and related transformation have immediate applications. 
With the aid of this method, the outstanding problem of the stability in phage lambda 
genetic switch has been solved Q]. More specifically, the relative stability of two fixed 
points in a dynamical system, in the presence of noise, can be uniquely determined. The 
connections of this mathematical structure to physics, and detailed behaviors near a fixed 
point, stable or unstable, have been analyzed rather exhaustively lfj. However, it is also 
noticed that in more general situations, the new transformation may give a result different 
from the conventional ones from the traditional i nteg ration of a SDE, such as by the Ito 
integration or by the Stratonovich integration <], The unique structure suggests yet 
another stochastic integration. The purpose of the present paper is to explore this suggestion 
and to put it into a statistical physics perspective. In this way the application of statistical 
physics to complex dynamics is explicitly suggested. 

In the following, we first summarize in section II the transformation from the standard 
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SDE into the transformed stochastic differential equation. We point out an interesting 
invariant of the transformation in section III, suggesting that the transformation might have 
a deeper mathematical root. In section IV we explicitly demonstrate the transformation in 
a special one dimensional case, where analytical results have been known in literature. A 
comparison to the Ito and Stratonovich integrations is made. Implications to various realistic 
applications are discussed in section V, and we conclude in section VI. Some mathematical 
details are provided in the Appendix. 



II. REVIEW OF THE TRANSFORMATION 



Let us consider a system described by state vector, q T = (qi,...,q n ), in n dimensional 
phase space, 1R™. Here the superscript r denotes the transpose. The dynamic equation 
describing the time evolution of the system is the usual Langevin like, standard SDE, 

q = f(q) + C(q,t) , (1) 

a shorthand writing for dqi/dt = fi(q) + d{q, t). Here ( T (q, t) = (Ci(<?, t), ■ ■ ■ , Cn(q, t)) denotes 
stochastic forces, assumed to be represented by Gaussian white noises, with zero mean, 
(C{q,t)) = 0, and variance 

(aq,t)C(q,t')) = 2D(q)5(t-t') , (2) 

where D is a positive semi-definite diffusion matrix. The angular brackets denote the average 
over noise distribution, and 5{t) is the Dirac delta function. The temperature in the present 
paper, whenever it can be defined, is always set to 1. This is equivalent to absorbing the 
temperature into the potential function to be discussed below. 

^, [h], [ll] that this equation may be transformed into the form 

[S(q) +r(g)]g = -V0(g) +£(</,*), (3) 



It was noticed 



where the matrix S is positive semi-definite and symmetric, the matrix T antisymmetric. 
In analogous to the overdamped Brownian system in physics, here — Sq plays the role of a 
dissipative force with a non- negative energy dissipation rate q T Sq/2. The matrix T plays a 
role similar to that of a magnetic field. The potential function is <f)(q). The new stochastic 
force £ T (g, t) = ...,£ n (q,t)), which has the same origin as i), is chosen to satisfy 

(aq,t)C(q,t')) = 2S(q)5(t-t') J (4) 



with zero mean. The steady state distribution for Eq. if exists, would be Boltzmann- 
Gibbs like 



The transformation from Eqs. (1) to (3) is rather remarkable in that the original SDE 
in Eq. (OQ), with no explicit potential function, is transformed into the dynamic equation 
in Eq. ([3]), governed by the potential <f), which can be obtained without actually solving 
the time-dependent equation. The usual assumption of zero probability current or detailed 
balance as the potential condition is not required in the passage from Eq. (pQ) to Eq. ([3]). 
As pointed out, the results obtained in this way rely on the partial differntial equation for 
the probability density, p(q,t), that is different in general from those based on either Ito or 
Stratonovich stochastic integration. Thus several questions naturally arise: Can this differ- 
ence be more clearly demonstrated? Are there further implications of this transformation? 
In the following we attempt to give an affirmative demonstration to those questions from 
two different angles: the invariant of the transformation and the treatments of SDE's. 

III. INVARIANT OF THE TRANSFORMATION 

Here we wish to point out an invariant for potential function during the transformation 
from Eq. (1) into (3) within the framework of Langevin equation, Eq.(6), and the Fokker- 
Planck equation (Klein-Kramers equation), Eq.(7), two of most important equations in the 
description of nonequilibrium processes in physics. This invariant suggests a generic nature 
of the present transformation deeply rooted in statistical physics. 

Let us consider the following nonlinear stochastic dynamics familiar in physics, an exam- 
ple for noise acting only on half of state variables: 



Here q, p are n-dimensional vectors, m a parameter, 5aiifixn semi-positive definite sym- 
metric matrix, T an n x n antisymmetric matrix, and <p a potential function. The subscript 
q in gradient indicates that it operates on q only. There is an additional relationship, due 
to Einstein, between the noise ( and S: (((q, t)( T (q, t')) = 2S(q)5(t — t'). This equation is 
the Kramers' dynamics in the standard Langevin form in physics [2|. 




(5) 




(6) 
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The probability distribution function p(q, t) corresponding to the dynamical equation 



121 ]. or Klein-Kramers equation [2|, 



satisfies the Fokker-Planck equation 

{d t + %-d q + J(q,p)-d p 

-d;S{q) [E. + d p }}p(q,p,t) = 0. 

Here f(q,p) = p T T(q)/m - <p q (q) with <p q (q) = d q <p(q). 



(7) 



A stationary state, while dtp = 0, exists and can usually be reached after a long time of 
dynamical evolution. From the analogy to physics, the final distribution from this dynamics 
is given by the Maxwell-Boltzmann-Gibbs distribution: 



p(q, p,t = oo) oc exp < — 



p- 

2m 



(8) 



There is a clear separation of variables q,p in the final distribution function. This suggests 
that the mass of the particle is not an essential quantify to determine the distribution in 
the g-space, a fact well known in classical statistical physics. Hence the validity of Eq. (5) 
is implied. 

We note that Eq. (6) is in the form of Eq. (1). To be more suggestive, we rewrite Eq. (6) 

as 

x = T{x) + Z(x, t) 



(9) 



x 




pjm 

■[S{q)+T(q))p/m-V q cf>{q) 



Here 

is a 2n-dimensional vector, 

F{x) = 

and 

Z 

The 2n x 2n diffusion matrix T> is defined by 

(Z(x, t)Z T (x, £')> = 2V(x)8(t - t') 

and is 

T>(x,t) 





(10) 



(11) 



The important question is whether we can obtain the corresponding form of Eq. (3): 

[S(x)+T(x)]x = -V x Q(x)+E(x,t) (12) 

with the same prediction on the potential, that is, without any change in physics prediction. 
Here T is a 2n x 2n antisymmetric matrix and S a 2n x 2n positive semi-definite symmetric 
matrix. We also require 

(E(x, t)E T {x, t')) = 2S(x)8(t - t') , (13) 

with zero mean of S. 

The answer is affirmative, demonstrating a very suggestive invariant of the transforma- 
tion. We give the result here, which can be verified directly: 

nx) J ™ J S f °). (14) 

Here / is the nxn identity matrix. The corresponding expression for V + Q = (S + T)' 1 is 

G(*) = [ ° " J | , (15) 
and V is precisely given by Eq.(ll). The potential is 

2 

Hx) = h + m ' (16) 

also precisely the same as that before the transformation. Thus, once the form of Eqs. (3) 
and (6) are reached, further manipulation of equations according to the transformation does 
not change the physical prediction. 



IV. TREATMENTS OF STOCHASTIC DIFFERENTIAL EQUATIONS 

The kinetic momentum in Eq. (6) can be eliminated, resulting in the following equation: 

mq = -iS(q) + T(q)}q - V,0(g) + £(q, t) . (17) 

Apparently, one would be able to take a zero mass limit by setting the right hand side of 
Eq. (17) zero. Carrying out this procedure, we have, copying Eq. (3), 

[S(q) + T(q)}q = -V q (f>{q) + £(q, t) . (18) 
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This can be converted into a more suggestive form, the standard form of SDE, copying 
Eq. (1): 

q = f(q) + ((q,t). (19) 

Here f(q) = -[%) + T(q))- 1 V q (f>{q) and C(g,t) = [S(q) + T(q)]^ 1 ^(q, t). 

Now, in connection of the singular behavior of the zero mass limit, we are facing a choice 
on the proper way to carry out the integration of Eq. (19). Different stochastic integration 
method would give different final distribution, hence different physical prediction. Here we 
illustrate this situation by considering explicitly the one dimensional case, n — 1. Higher 
dimensional case will be considered elsewhere 29j. In one dimensional case, n = 1 and 
T = 0, and Eq. (3) or (19) becomes 

q=-D(q)V<Kq)+C(q,t), (20) 
and S(q)D(q) = 1. We explore its three different integrations below. 

I. Present process. If we follow a limiting process to reproduce the final steady state 
distribution, the dynamical equation in this limit would be 

d t p(q, t) = [d q D(q)d q + d q D{q)<P q (q)]p(q, t) . (21) 



This is in the form of the celebrated Smoluchowski equation [13j. Three different deriva- 
tions from Eq. (7) to Eq. (21) are provided in the Appendix. Thus, the present stochastic 
integration process overlaps with that of the Smoluchowski process of strong damping limit. 
The steady state distribution would be, if exists, 

p(q, t = oo) oc exp{-0(g)} , (22) 

the same as in Eq. (5). 

II. Ito process. The Ito stochastic integration [ijj is a strict implementation of Markov 
process and a martingale: There is no memory effect of previous dynamics and there is no 
information from future. The usual differentiation and integration rules would not be directly 
applicable here. Fortunately, this case has been studied in great details. The corresponding 



Fokker-Planck equation is, from Eq. (20) 15] : 



d tPl (q, t) = [d q d q D(q, t) + d q D(q, t)<t> q {q)\ Pl {q, t) . (23) 



The subscript / indicates that Eq. (23) is in accordance with Ito process. The steady state 
distribution would be, if exists, The steady state distribution would be, if exists, 

p 7 (g,t = oo) oc -^yexp{-0(g)} . (24) 



III. Stratonovich process. In the Stratonovich stochastic integration, the usual dif- 



ferentiation and integration rules apply 16j. This implies that the stochastic process such 



specified is a Markov process but not a martingale. Again, the stochastic integration has 



been carried out, and the corresponding Fokker-Planck equation is, from Eq. (20) 171 ]: 

d tPs {q,t) = [d q D*(q,t)d q Di{q,t) 

(25) 

+d q D(q,t)<P q (q)]p s (q,t) . 

Similarly, the subscript S indicates that Eq. (25) is in accordance with Stratonovich process. 
The steady state distribution would be, if exists, the steady state distribution would be, if 
exists, 

ps(q, t = oo) oc } exp{-<p(q)} . (26) 
Da(q) 

The above one-dimensional examples clearly show that the present treatment of SDE 
implied in the transformation may be regarded as another stochastic integration process. It 
can indeed produce results different from those of Ito or Stratonovich process. One way to 
understand this difference is to view the zero mass limit as a singular limit in which specific 
procedure must be identified in order to get a well-defined result. As explicitly demonstrated 
above, there are relative shifts of density peaks from one process to another. In addition, 
the positions of peaks of the distribution function may not coincide with those of zeros of 
the potential cf) or the fixed points of /. We have observed such shifts during numerical and 



analytical study even when the diffusion matrix D is constant 18J, that is, even when there 



is no difference among those of Ito and Stratonovitch, and more general considerations 
Similar shifts have also been observed before 



20 



19] 



2l| . They are associated with the absence of 



detailed balance: In the presence of detailed balance condition and with a constant diffusion 
matrix, all those stochastic treatments would lead to the same potential function. 



S 



V. DISCUSSIONS 



With various limiting procedures to integrate the SDE of the form of Eq. (1), (19) or 
(20), which one would be right? The answer is, not surprisingly, the procedure has to be 
determined by the real problem at hand and any of them can be correct in a practical 



situation 
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221 ]. We illustrate this point by examples. 



It should emphasized that it has been shown mathematically that in general, non-white 
noise can be constructed and its limiting process can lead to different stochastic integration 



with different Wong-Zakai correction 



191 ] . We also note, in the case of Id, that the Eq. fT22j) 



corresponds to the choice of evaluating the function at the end of a time interval. Higher 



dimensional case would be different 



291 ] . The Ito process, Eq.(23), is then a choice of initial 



point, without peering into future. The Stratonovich process, Eq.(25), is a choice of middle 
point. 

It is evident that from the physics point of view the first process appears to be a natural 
choice: the presence of small mass limit encountered frequently. It results in the known 

nnn 

Smoluchowski equation [2|, |7|, [12[ . Numerous experiments in physics exist to support this 
approach. We wish to point out that the form of Eq. (6) allows a generalization into 
more complicated situations, such as the colored noise, as done in the dissipative dynamics 



23. 



241 ] . This is also the form used in a recent development of a mesoscopic, open-system 



nonequilibrium thermodynamics 25 1. 

From an engineering point of view, everything would be thought as a limit of continuous 
process, hence Stratonovich process would be a natural choice. Indeed, analog experiments 
have been done to check this situation. The results are in remarkable agreement with the 
Stratonovich stochastic integration 26]. 



From a population geneticist's point of view, Ito process would appear more natural to 



27|, 



28|. We 



model the population dynamics. Indeed this view has been carried out in detail 
are, however, not aware of any precise comparison between experimental/expirical data and 
the theoretical calculation which can single out the Ito process as the appropriate description 
in this case. 

It is worthwhile to point out here that, in connection with the invariant discussed in 
section III, in the corresponding Fokker-Planck equation for Eq. (6), the Eq. (7), there is 
no need to differentiate all those stochastic processes discussed in section IV, because of the 
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special algebraic structure expressed in Eq. (11) and (14). This adds a further preference to 
view the present process as a "natural" one. 



VI. CONCLUSION 



In this paper we have given an explicit demonstration of the difference between the present 
integration of the stochastic differential equation and those according to classical methods in 
one dimension. While we believe all those integration methods are in themselves consistent, 
we do point point that the present process has a certain advantage, as demonstration by the 
invariant of the transformation and by the zero mass limit. Thus, the existence of potential 
landscape in complex systems is put on a firmer theoretical ground. 

Acknowledgments. Discussions with R. Mannella, P. McClintock, D.J. Thouless, D.B. 
Walton, L. Yin, and X.-M. Zhu are highly appreciated. This work was supported in part by 
the USA National Institutes of Health through Grant number HG002894 (PA.). 



Appendix: Three Derivations of Eq. (21) 

The particular form of the partial differential equation in Eq.(21) has not been widely 
discussed in literature. Since it plays an important role in our study of the potential land- 
scape, we give three derivations which are not related to the usual treatments of stochastic 
differential equations 19|. An extension of Eq.(21) into higher dimensions is discussed in 
Ref.LoL 



First derivation 

Here we give some justifications, both physical and mathematical, to the derivation of 
the Eq. (21) in one dimension. The derivations for the equations (23) and (25) are well 
documented in the literature. Interested readers should consult textbooks such as that by 
van Kampen jjj, Gardiner 7], or Risken [l^ . 

We note that there is an alternative expression for Eq. ([6l) 

Q = v 

(27) 

mv = -[S(q,t)+T{q,t)]v-V q <Kq)+({q,t) 
Eqs. (jSJ) and ff27l) are for the inner and outer solutions to the problem with fast and slow 
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time scales respectively. In one dimensional case, n = 1 and T = 0, we may rewrite equation 
(7) as 

9 t p(q,P,t) = {-%dg + (pg{q)d p 

+ d p S(q) [E.+d p ])p(q,p,t). 
The right hand side of Eq. (28) can be regrouped according to the following manner: 



(28) 



(29) 



d t p{q,p,t) 

= ( dp [S(q)d p + S(q)i\+d p( p q (q)-i9 q ) 

= (fc- a *35j] ism+s(q)n 

+d,sk [S(q)d p + %)£] + d P Mq) - %d q ) p(q,p,t) 

= ([ d p- d «W)] [S(Q)d P + S(q)^\ 
+d p [d q + 4> g (q)}) p(q,p,t) 

= ([ 9 p- d «m] l s ^ 9 p + s ^ + + 9 «] 

There are two lines of attack for this problem in the literature. One was taken by Kramers 
himself and Chandrasekhar [1], the other is based on the projection operator approach Q. 
The essential assumption in the latter approach is that the momentum p stationarity is 
achieved instantaneously when the mass m is small. It takes the time of order m/S(q) 
for it to reach equilibrium. Because we are interested in the dynamics on the time scale 
At ^> m/S(q) , we may regard the dynamics of the kinetic momentum p is in an instantaneous 
equilibrium with the dynamics of the coordinate q. We may then assume the trial solution 
for the distribution function in the following form: 



p(q,P,t) = p(q,t)p K (p,t\q) 



(30) 



with conditional probability of p given q: 
pK(p,t\q) = 



2ixm 



cxp 



2m 



P 



S(q) 



(31) 



Here A(fi(q) = — lnp(g, t) — (f>(q) and A(j) q (q) = d q A(j)(q). The distribution of kinetic mo- 
mentum is Gaussian, in the familiar form of kinetic energy. At a given moment the average 
kinetic momentum is not zero, corresponding to the drifting velocity A<fi q (q) / S(q). The 
drifting velocity is zero at the eventual equilibrium state, as would expected. We expect 
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that the error raised from this choice will be in the higher order of m/S(q) comparing to 
what would be kept. We will show below that it is indeed the case. 

With the trial solution the contribution of the first term of the right hand side of Eq. (29) 

is: 



K 



d - d — 

u p u i S(q) 



x 



X 



[S(q)d p + S(q)% + (f> q (q) + d q ] p(q, t)p K (p, t\q) 
p(q,t) \d p - (-A^(g) - <t> q {q) + d q )^ 
[S(q)d p + S(q)^ + ( f )q (q) 
+(-A<f> q (q) - (f> q (q) + d q )} p K (p, t\q) 
p(q,t) d p - (~A(f) q (q) -4> q (q) + d q )^ 
[S(q)d p + S(q)% - A<f> q (q) + d q )] p K (p, t\q) 



= P(q,t) [d P - {-A<j)g{q) - (f) q (q) + d q )^ 
[S^dp + S^-A^M + d^x 

2 



X 



1 



2-m eX P 1 2m [P S(q) 

Moving pk from the right to left, and grouping all terms, we have 



mAtj> q (q) 



K = P(<l,t) x y^expj-^ (p- 1^1) 2 } 
"i(f-1f)-(-A^)-^) 

mA(f> q (q) 



X 



d„ 



■ mA<t> q (q) 



-Sfe)i(p- 
-i(p- 



S(q) J "</V S(q) 
mA<t> q (q) 



+ 9q))sh 



X 



S{q)i-A<P q {q) 



mA4> q (q) 
S(q) 



n 



■ rnA<t> q (q) ' 



S(q) 



P 



mA4> q (q) 
S(q) 

mA0q(g) < 



= pM x \/S exp 

_1_ ( n _ mA<j> q {q) \ ( ', , i p. /mj s . :,,':/ , 
m [P S(q) J \ \ Q\ S(q) >> 

-d q )sk)) - (-A^(g) - U<l)) 

J_ („ _ mAcl> q (q) \ a , mA<j> q \ 
"m \P S(q) ) l\ S (q) 



X 



(32) 



(33) 
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Performing integration over p, all linear terms are zero, we have 




(34) 



when — > 0. Note that p(q,t) = J^dp p(q,p,t), integrating Eq. (29) over p, we have 



This is Eq. (21). The correction is in higher order of m/S. This completes our first derivation 



Second derivation 

We can also follow the original approach of Kramers for the derivation of Eq. (135]) . our 
second derivation of Eq.(21). One subtle issue in this approach is realizing that to obtain 
p(q,t) from p(q,p,t), one needs to integrate out p along the curve q + p/S = constant 
in which (p/m) is the velocity and m/S is the relaxation time for the velocity. For any 
v on the order 1, of course, this correction is on the order of m/S. For non-constant 
S(q), the derivation presented below is in close parallel to that of Kramers, but with extra 
mathematical complications. 

We first observe several mathematical relations. 
(A.i) From distribution p(x, y), the distribution for z = x + ay is J p(z — ay, y)dy. 
(A.ii) We have an identity 




(35) 



of Eq.(21). 




(36) 



d^(q 



' m 



z 



Til 



1 



JS(q)dq) 



in which ^(g, v ) 
(A.iii) Let v — - 



rn 



= —m^/(q,v)/S(x). 

(z — J S(q)dqj and its inverse function be q = g{z — mv). Then we have, 
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for any G 



l r°° 

m 



S^d q G( q ^(z-jS{ q )d q )) 

-oo , G{g(z-mv),v) 
oo S(g(z—mv)) 

9=9(2) 



(37) 



f oo 
-oo 



" U 5(g) 



when m — > 0. 

Starting with the SDE in Eq. (1271) . we have equation similar to (29): 



d t p(q,v,t) 



d - d m 



S(q) 



d v + 



vS(q)+<j> q (q) , 1 



1 S(g) _ 

+ d g ^[<9 g + 9 (g)])p(g,M) • 



9 a 

m y 



(38) 



Define random variable z = mv + f S(q)dq and we have its distribution according to the 
relation (A.i), 



p{z,t) 



m 



z 1 

P 1 9, 

m m 



S(q)dq dq. 



(39) 



Using the above relations (A.ii) and (A.iii) 

d t p{z,t) 

1 c +oo , 

= — ctox 

m J — oo ^ 



I^^[d q {sbd q + ^) P (q,v,t) 



9=9(2) 



(40) 



1 

S(q) 



q=g(z) 



Finally, we note that if we define random variable ( = g(u), then the distribution for (: 
P((y) = S(y)p z (J S(y)dy). In the limit of m — > 0, ( = g (mv + J S(q)dqj — > q, and the 
distribution for ( satisfies the equation: 



d t p(q, t) = d q 



-dq + 



(?) 



S(q) « S(q) 



p(q,t). 



(41) 



Again, this is Eq. (21). 



Third derivation 

It may appear that the first derivation is too intuitive and the second one is too formal. 
One may wonder whether or not there is a more programmatic derivation. The answer is 
positive. Here we give such one where the higher order correction in the small mass limit 
would be systematically obtained. 
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The trial distribution function is assumed in Eqs. (30), (31). We show it is indeed the 
right solution from the Klein-Kramers equation in Eq. (28). 
First we write 

,2 



V 
2m 



p(q,p,t) = exp ( - A(j){q) - 

+A(q,p,t) + mB(q,p,t) 

The higher order terms in m are included in B(q,p,t). Then Eq. (28) gives 



(42) 



p_ 

m 



l h + At + mB t 
-(A(j> q - A q ) - pB q + (j) q A p + m<f) q B p + S(q) 

(-—A p - pB p + A 2 p + 2mA p B p + m 2 B^ (43) 

where the subscripts denote the derivatives with respect to t, q,p. Then the term of 0(l/m) 
yields 

— A q — S(q)A p = . (44) 
Then we can find A for small p, noting p is order of m. 



A(q,p) 



(45) 



This satisfies Eq. (1441) in the above, where A q gives 0(p). This justifies the trial form of 
the distribution function used in our first derivation. The rigorous solution for A(q,p) can 
be found: 



A(q,p) 



1 — exp 



P a 
S{q) dq 



A0 



(46) 



We will take only the first order in p. Nevertheless, we point out that this procedure can be 
carried out to higher orders by matching the corresponding orders in Eq. (43). 



15 



The expression corresponding to Eq. (33) should read as follows: 



K = p(q,t) x 



cxp 

2ixm 2m 



1 / _ mA(f) q (q) 



x 



-— ( /; 

m 



+ —[P- 

m 



mA(p q (q) 

S(q) 
mA<p q (q) 



S(q) 
■ &<t> q {q) ~ <f> q (q) 



mA<p q (q) 
S(q) 



}s(q) 



x 



m \ S(q) J m 



-A6 (a) + 1 ( p _ m *MQ) \ d ( mAMg) 



p(q,t) x 



1 



exp 

2nm 2m 



I ( mA(j) q (q) 



x 



I ( mA(j) q (q) 



m 



1 



S{q) 
s 1 



S(q) 



S(qr q \S(q) 



x < — I p 

m 



qJ S(q) 
mA(j) q (q) 



S(q) 



\S(q)J 



According to this, Eq. (34) leads to: 



d, 



mA(j) q 



S(q) 

+ q \S(q)J 
. 



+ 



m_ f / AtA 
S(q) V q \S(q)) 



(47) 



(48) 



It is exactly zero, not approximate zero as in Eq. (34) Again, Eq. (21) follows in the zero 
mass limit. 

Above demonstration evidently indicates that the novel stochastic integration procedure 
discussed in the present paper has a strong connection to the overdamping dynamics, the 
Smoluchowski limit. They are the same in one dimension, but differ generally in higher 
dimensional cases which would be addressed elsewhere. 
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